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ABSTRACT 



We investigated the dependence of the solar magnetic parity between the 
hemispheres on two important parameters, the turbulent diffusivity and the 
meridional flow, by means of axisymmetric kinematic dynamo simulations based 
on the flux-transport dynamo model. It is known that the coupling of the mag- 
netic field between hemispheres due to turbulent diffusivity is an important factor 
for the solar parity issue, but the detailed criterion for the generation of the dipole 
field has not been investigated. Our conclusions are as follows. (1) The stronger 
diffusivity near the surface is more likely to cause the magnetic field to be a 
dipole. (2) The thinner layer of the strong diffusivity near the surface is also 
more apt to generate a dipolar magnetic field. (3) The faster meridional flow is 
more prone to cause the magnetic field to be a quadrupole, i.e., symmetric about 
the equator. These results show that turbulent diffusivity and meridional flow 
are crucial for the configuration of the solar global magnetic field. 

Subject headings: Sun: activity — Sun: interior — Sun: dynamo 
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Introduction 



The solar magnetic 



internal ionized plasma ([Parker 



1-year 



cycle is thought to be sustained by the dynamo motion of the 



19551 ). Based on the internal structure of the velocity field, 



i.e. , the meridional 



iow a nd the differential rotation revealed by heli oseismology (see review 



Thompson et al. 



Dikpati fc Charbonneau 



2003 ), the flux-transport dynamo was suggested ( Choudhuri et al. 



1999 



Kiiker et al. 



2001 



Hotta fc Yokoyama 



1995 



20101), as a model to 



successfully explain some features of solar activity such as the equatorward migration of 
sunspots and the poleward migration of the surface field. 

The solar global field has a distinct parity: the poloidal field is a dipole, i.e., 
antisymmetric about the equator. The polar fields almost always have the different sign 
between hemispheres, even though they show the occasional weak north-south asymmetry 
in phase and amplitude. In addition, Hale's polarity law states t hat the sun spots between 



19081 ). It can then 



hemispheres are nearly always antisymmetric about the equator (iHald 
be interpreted that the toroidal fields (-B,^) below the surface are antisymmetric about 
the equator. This interesting feature is, however, not axiomatically explained by the flux 
transport dynamo model since this model significantly depends on three free parameters, 
i.e., the a-effect, the meridional flow, and the turbulent diffusivity. 



It has been suggested that the a-effect around 



;he base of the convection zone leads to 



the generation of the g 



2002 



Chatter] ee et al. 



obal dipolar magnetic field. (jPikpati fc Gilman 



2001 



Bonanno et al. 



20041 ). The existence of the poloidal fields around the tachocline and 



the coupling of these fields between hemispheres are significant factors for the generation 



of the dipole fie 



Chatterjee et al. 



d . A d etailed explanation of this process is given in the next paragraph. 



(120041 ) also suggested, however, that the dipole field can be obtained with 
the strong diffusivity in the convection zone without the presence of the the a-effect around 
the tachocline. Hence the exact necessity of the a-effect in generating the dipole field is still 
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inconclusive. 

The dependence of the parity on these parameters can be explained when we understand 
the role of the turbulent diffusivity in the solar magnetic parity issue. If the global magnetic 
field is antisymmetric, i.e. is a dipole like our sun, the component of the magnetic vector 
potential in each hemisphere has the same sign (Fig. [T^). When the cyclic phase in one 
hemisphere slightly differs from the other, the coupling effect by the turbulent diffusivity of 
the poloidal field distinguishes the phase difference in the vector potential and causes the 
magnetic field to be a dipole. On the other hand, when the magnetic field is symmetric, 
i.e., is a quadrupole, this effect does not occur. Therefore the substantial coupling of 
the poloidal field generates the antisymmetric (a dipole) magnetic field. The sign of the 
toroidal field in one hemisphere differs from that in the other hemisphere when the global 
magnetic field is a solar-like dipole. With the same logic posited above, it is obvious that 
the substantial diffusive coupling of the toroidal field between the hemispheres helps the 
magnetic field to become symmetric (a quadrupole; Fig[T|D). In summary, the parity of 
the stellar global magnetic field depends on which field, the toroidal or the poloidal, is 
more coupled by the turbulent diffusivity between the hemispheres. Detailed systematic 
parameter studies are necessary to understand for the parity issue. 

In this study, we investigate the dependence of the solar magnetic parity on the 
distribution of the turbulent diffusivity and the amplitude of the meridional flow. The 
obtained constraint on the turbulent diffusivity is important since it is one of the key 
components of the solar dynamo model, although it is difficult to measure by direct 
observations. 
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Model 



We solve axisymmetric kinematic dynamo equations. The magnetic field is divided into 
the toroidal field B and the poloidal field Bp = V x [A{r, 9)e^] in the spherical coordinate 

as 



B = B{r,9)e^ 



Br 



where A is the component of the magnetic vector potential and is the unit vector along 
the direction. Then the following standard forms of the dynamo equation are derived 
from the magnetic induction equation as 



dB 1 
dt r 



g^irUrB) + -{ueB) 



rsin6'(Bp ■ V)Vt - {Vrj x V x Be^ 



(2) 



dA 1 
dt r sin 9 



1 



5, 



sin^ 9 ^ 
(u ■ V)(rsin6'A) 



A + S{r,9;B) 



(3) 



sin^ 9 ^ 

We specify the meridional fiow u = w^Cr + ueeg, the differential rotation Q, and the 
turbulent diffusivity rj. A source term S{r, 9; B) is artificially added to the right-hand 
side of equation ([3]). This source term describes the generation of poloidal fields at the 



(Leighton 


1969; 


Bab cock 


1961) 



196ll ). Once these quantities are specified, we can solve equations 



(|2]) and ([3]) to study the behavior of the dynamo kinematically. 



The formula for the differential rotation is given as 



n{r,9) = n^ + 



1 + erf 2 



dr 



{n,{9) - fij 



(4) 



where VLs{9) = Q-Eq + a2 cos^ 6' + 04 cos^ 9 is the surface latitudinal differential rotation and erf 
is the error function. The parameters are set as ^c/Stt = 432.8 nHz, Q-Eq/2n = 460.7 nHz, 
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a2/27r = —62.69 nHz, aA^jln = —67.13 nHz, and Vc = 0.7R and dc = 0.05i?, which closely 
resemble the best-fit helioseismic solution. R is the solar radius, and is the central 
radius and the thickness of the tachocline, respectively. Qc is the rotating rate of the 
core. This differential rotation profile has a purely latitudinal difference with equatorial 
acceleration in the convection zone. It smoothly matches across the "tachocline" with the 
core rotating rigidly. 



We now describe how the meridional flow is specified. While a poleward meridional 
flow is observed near the solar surface, the structure of the internal return flow is at 



present unconstrained observa t ional. 



y. We choose here an analytical form suggested by 



van Ballegooijen fc Choudhuril (119881 ) with the density proflle in the convection zone given 



by p(r) oc [{R/r) — 1]"^; the r and 9 components of this flow are as follows: 



where 



X 



1 



+ 



Cl 



C2 



m + 1 ' 2m + 1 2m + p + 1 
X esin''^[(g + 2) cos^^-sin^^] 

3 



m+p 



ue(r, 



Uo fR 
" f \r 

X sin'^'*'"'^ 6 cos 6 



1 + cir - 



m+pl 



Cl 



ar) = - - 1, 

r 

(2m + l)(m + p) 
(m + l)p 



C2 



(2m + p + l)m^_(^+p) 



(m + l)p 



(5) 



(6) 



(7) 



(9) 



1. 



(10) 
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As shown in eqs. ([5]) and Qj, uq is the amphtude of the velocity and p and q are respectively 
the radial and latitudinal dependence of the flow. We specify the bottom of the meridional 
flow Tp, and the normalization constant / to set the maximum speed of the meridional flow 
of the 6 component to uq- We use the parameter values m = 0.5, p = 0.25, g = and 

flow slightly penetrates into the rigidly rotating core. 



Vp = 0.62R. Note that our mer idiona. 



i.e. Tp < Tc (iHotta fc Yokoyama 



2010h . 



We assume that the net magnetic diffusivity in the convection zone is dominated by its 
turbulent contribution. We adopt a diffusivity profile of the form 



=Vcore + 



Vt 



1 + erf 



+ 



1 + erf 



r — r2 
do 



(11) 



Here, ri = 0.7i?, di = 0.02R and d2 = 0.02R. This profile consists of three layers. In the 
strong diffus i vity l ayer (r > r2), the diffusivity is prescribed by rj^ within 10^^ — 10^^ cm^ s~^. 



Wang et al. 



argued that the surface diffusivity should be 6 x 10 cm s to be 



consistent with observations of the time development of the surface magnetic field. In the 
convection zone we use the fixed value rji = 5 x 10^° cm^ s^^. In the subadiabatically 
stratified core there is no turbulence (or at least, far less), so that the diffusivity is 
presumably much weaker. We use the value r/core = 5 x 10® cm^ s^^. For convenience, we 
define the surface depth dg = R — r2 which denotes the thickness of the strong diffusivity 
layer. We take dg and ?7s as free parameters, since they affect the parity of the magnetic 
field. 

Some parts of the toroidal field in the tachocline rise to the surface due to the magnetic 
buoyancy and generate active regions. During this process, the flux tube ex pands and 



the C oriolis force bends the flux. There is observational evidence offered by 



Babcock 



(119591 ) that the decay of tilted bipolar acti ye regions can produce a substantial amount 



of the net poloidal fields near the surface (jWang fc Sheeley 



199l|). This is called the 
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"Babcock-Leighton a-effect." Based on this, we assume that the poloidal source term is 
taken in the form 



X 



l + {B{r,,e)/B,,y 



'12) 



where 



a{r,e) =1 



erf 



r — r4 



erf 



r — rs 



(13) 



X sin 9 cos 9 



1 



1 _^ e"'^(^"'"/'^) 

The parameters are r4 = 0.95/2, = R, d4 = 0.05/2, d^ = 0.01/2, and 7 = 30. We 
concentrate the a-effect by the last factor in eg . ([I3ll at the low latitude in which there are 



many observed active regions ( jPikpati et al. 



20041 ). The source term is made proportional 
to the toroidal field strength at the same latitude in the tachocline r = Tc (eq. dl])), i.e. the 
base of the convection zone, since it is assumed here to originate from the radially emerged 
magnetic fluxes. The quenching term {1 + [/?(rc, ^^)//?cq]^}^^ in eq. f|T2|) ensures that the 



poloidal fleld production r apidly vanishes as the deep toroidal fleld strength exceeds B, 



( ICattaneo &: Hughes 



eq 



19961 ). /?eq is the equipartition magnetic fleld. The Coriolis force 
cannot bend a strong magnetic fleld (> /?eq) in the convection zone. We use a flxed value 
/?eq = 4 X 10*^ G. So is determined by flxing the value of the dynamo number (see ^ in 
detail) . 

We solve equations (|2]) and (|3]) numerically in all the sphere of the meridional plane 
in 0.6/2 < r < R and < < vr with the modifled Lax-Wendroff scheme. We use a 
moderate resolution of around 64 grid points in the radial direction and 128 grid points 
in the latitudinal direction. At the lower boundary (r = 0.6/2), we set both B and A 
at zero, indicating that the radiative core is a perfect conductor. At the top boundary 
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(r = R), we set -5 = and smoothly match A onto an exterior potential field solution 
( Dikpati fc Choudhuri Il994 ). At both poles {6 = and tt), we set B = A = for the 



regularity. The numerical convergence is checked by runs with different grid spacings. 



3. Results 

A new indicator of the magnetic parity is defined in this study. The radial magnetic 
field at the surface can be decomposed as 



Br{R,e)=J2(^nPn{cOSe), (14) 



n=0 

where P„ is the Legendre polynomial. Then we define the symmetric parameter as 



1=0 i=0 



C2i+1 1 



SP = ^ . (15) 



E 



Ci 



i=0 

Each even (odd) order of the Legendre polynomial is symmetric (antisymmetric) about the 
equator. Therefore, SP = 1 corresponds to the purely symmetric mode about the equator 
and SP = — 1 is the antisymmetric mode. 

We first show a representative reference solution, computed with the amplitude of 
the meridional fiow uq = 1000 cm s~^, the amplitude of the a-effect sq = 160 cm s~^, 
the amplitude of the turbulent diffusivity near the surface layer rj^ = 2 x 10^^ cm^ s~^ 
and the thickness of the strong diffusivity layer dg = O.IR. The results are shown in the 
time-latitude plots in Fig. [2l This simulation is started with a symmetric initial condition. 
As time passes, the global magnetic field becomes antisymmetric. In such a case, the 
symmetric parameter develops as shown in Fig. [31 The black (red) line denotes the result 
with a symmetric (antisymmetric) initial condition. Regardless of the initial condition, 
the parity of the magnetic field approaches the stationary antisymmetric state where the 
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symmetric parameter becomes ~ —1. We conclude that the magnetic field with these 
reference parameters finally becomes a dipole field. 

We investigate the asymptotic stationary values of the symmetric parameter for runs 
in different setups. We carried out runs by choosing a value for the surface diffusivity ?7s, 
from 8 points in the range 6 x 10^^ to 1 x 10^^ cm^ s~^ and the surface depth dg, from 5 
points in the range 0.1-R to 0.25-R. We specify the amplitude of the a-effect by 

So = 160 cm s"^ ( t!^ — ^ (-^^ (16) 

V2 X 1012 cm2 s-iy \0.1RJ ^ ' 

The background reason for this setup is to fix the value of the surface integrated dynamo 
number, 

where and k denote the wavenumber of the magnetic field and da denotes the thickness of 
the layer where the a-effect is effective, and where k^-, k and d^ are assumed to be unchanged 
from case to case. This dynamo number is kept unchanged since it determines the oscillatory 
nature of the dynamo. The reason for this idea is because the strong a-effect is necessary 
to endure the diffusivity in the large area. We also vary the amplitude of the meridional 
flow: mq = 1000 cm s~i(slow meridional flow case) and 2000 cm s~i(fast meridional flow 
case). For every parameter set, we conducted runs with both symmetric and antisymmetric 
initial conditions to ensure that the asymptotic value of the symmetric parameter does not 
depend on the initial parity. 8 (diffusivity) x 4 (surface depth) x 2 (meridional flow) x 2 (initial 
parity) = 128 simulation runs carried out. All simulations are calculated for more than 
10000 years. 

There are two types of solutions when the value of the SP is around zero and we 
categorize them as "mixed-partity" cases. One type is similar to the reference case (Fig. 
UK). The value of the SP flnally converges. The other type is interesting in that the value 
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of the symmetric parameter does not converge and continues to oscillate between the 
quadrupole and the dipole solutions (Fig. HJo). Since the averaged value in the calculation 
duration is close to zero, we adopted it for the SP in such cases. 

The results of this parameter space study are shown in Fig. [51 The dynamo cycle 
period is also shown by the contour lines. The period is shorter when the surface depth is 
thicker since the transport of the magnetic flux by the diffusivity is more effective. Panel 
(a) shows the result of the slow meridional flow case [uq = 1000 cm s^^). Two points can 
be ascertained from this figure. One is that regardless of the surface depth, the strong 
diffusivity (> 3 x 10^^ cm^ s~^) can make the magnetic field to become a dipole (SP~ — 1). 
The other is that with the thinner surface depth, no strong diffusivity (> 1 x 10^^ cm^ s^^) 
is needed to generate the dipole field. This means that the magnetic field is more likely to 
be a dipole with the thinner surface depth. Fig. 4b shows the result of the fast meridional 
flow case {uq = 2000 cm s~^). It is obvious that the parameter area for the symmetric 
solutions, i.e. Sp > 0, increases. This indicates that the fast meridional flow causes the 
magnetic field to be symmetric. 

4. Discussion and Conclusion 

We investigated the dependence of the global magnetic parity on the distribution of the 
diffusivity (the amplitude and the surface depth) and the amplitude of the meridional flow. 
Three results were obtained. First, the model shows that the stronger diffusivity near the 
surface acts to make the magnetic field a dipole. The diffusivity near the surface enhances 
mainly the coupling of the poloidal field near the surface between the hemispheres, leading 
to the generation of dipolar magnetic field. The second result is that the thinner layer of 
the strong surface diffusivity also works to cause the magnetic field to become dipolar. The 
thinner surface depth suppresses the coupling of the toroidal field between the hemispheres 
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since most of the toroidal field exists around the tachocline. The third result is that the 
fast meridional flow causes the magnetic field to become a quadrupole. The fast meridional 
flow prevents the poloidal field from coupling near the surface of the equator because the 
flow transports the poloidal field poleward. In addition, the flow transports the toroidal 
field around the tachocline equatorward, thus causing the coupling of the toroidal field. 
These three results quantitatively constrain the distribution and the amplitude of turbulent 
diffusivity, which cannot be determined by observation and is a important factor for the 
dynamo problem. 

In this study, we did not investigate the dependence of the parity on the a-effect 
in the convection zone, which may be a strong factor in causing the magnetic field to 
become a dipole. The poloidal field generated by this effect around the tachocline is 
transported equatorward by the meridional flow, and this process enhances the coupling o f 



the poloidal field betw een the hemispheres (jPikpati fc Gilman 



Chatter] ee et al. 



2001 



Bonanno et al. 



2002 



20041 ). It is possible that the criterion for a dipole field we obtain in this 



study may be modified with this type of a-effect. We will discuss the possibility of the 
existence and the infiuence of the a-effect in a forthcoming paper. Another interesting issue 
to be addressed is the possibility that the variation of the velocity field in the solar cycle 



affects the parity. In the calculations for the earth dynamo the re is the significant di: 



between the kinematic and the MHD cases in the parity issue (INishikawa fc Kusano 



erence 



2008h. 



Thus , in the future we will investigate the parity issue with the Lorentz feedback (IRempel 
20061). 



Numerical computations were carried out on the general-purpose PC farm at the 
Center for Computational Astrophysics (CfCA) of the National Astronomical Observatory 
of Japan. The page charge for this paper is supported by CfCA. We have greatly benefited 
from the proofreading/editing assistance from the GCOE program. 



- 13 - 
REFERENCES 

Babcock, H. D. 1959, ApJ, 130, 364 
Babcock, H. W. 1961, ApJ, 133, 572 

Bonanno, A., Elstner, D., RMiger, G., & Belvedere, G. 2002, A&A, 390, 673 
Cattaneo, F., & Hughes, D. W. 1996, Phys. Rev. E, 54, R4532 
Chatterjee, P., Nandy, D., & Choudhuri, A. R. 2004, A&A, 427, 1019 
Choudhuri, A. R., Schussler, M., & Dikpati, M. 1995, A&A, 303, L29+ 
Dikpati, M., & Charbonneau, P. 1999, ApJ, 518, 508 
Dikpati, M., & Choudhuri, A. R. 1994, A&A, 291, 975 

Dikpati, M., de Toma, G., Gilman, P. A., Arge, C. N., & White, O. R. 2004, ApJ, 601, 1136 

Dikpati, M., & Gilman, P. A. 2001, ApJ, 559, 428 

Hale, G. E. 1908, ApJ, 28, 315 

Hotta, H., & Yokoyama, T. 2010, ApJ, 709, 1009 

Kiiker, M., Riidiger, G., & Schultz, M. 2001, A&A, 374, 301 

Leighton, R. B. 1969, ApJ, 156, 1 

Nishikawa, N., & Kusano, K. 2008, Physics of Plasmas, 15, 082903 
Parker, E. N. 1955, ApJ, 122, 293 
Rempel, M. 2006, ApJ, 647, 662 



- 14 - 

Thompson, M. J., Christensen-Dalsgaard, J., Miesch, M. S., & Toomre, J. 2003, ARA&A, 
41, 599 

van Ballegooijen, A. A., & Choudhuri, A. R. 1988, ApJ, 333, 965 
Wang, Y., Nash, A. G., & Sheeley, Jr., N. R. 1989, ApJ, 347, 529 
Wang, Y., & Sheeley, Jr., N. R. 1991, ApJ, 375, 761 



This manuscript was prepared with the A AS macros v5.2. 



- 15 - 




Fig. 1. — Illustration of the parity issue. Panel (a) shows the poloidal fields (line) for a 
dipole and a quadrupole field and the corresponding vector potentials. Panel (b) shows the 
toroidal field for a dipole and a quadrupole. 
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Fig. 2. — Butterfly diagram for tlie reference solution. Time-latitude plot of B^\r=o. 
contour is superposed on the color map of the surface radial fields. 
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Fig. 3. — Time-development of the symmetric parameters. The black (red) line corresponds 
to the results of the symmetric (antisymmetric) initial condition. Regardless of the initial 
condition the symmetric parameter finally becomes ~ — 1 (antisymmetric solution). 
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(a) convergent type 




(b) oscillatory type 




Fig. 4. — Two types in "zero" symmetric parameter cases. The time- development of the 
symmetric parameter for each type is shown. The format is the same as in Fig. [31 Panel (a) 
shows the convergent type in which the value of the symmetric parameter finally converges 
to zero. Panel (b) shows the oscillation type in which the symmetric parameter continues to 
oscillate between the quadrupole (SP~ —1) and the dipole (SP~ 1) solutions. 
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0.6 1.0 2.0 6.0 10.0 0.6 1.0 2.0 6.0 10.0 

diffusivity (10^^ cm^ s"') diffusivity (10^^ cm^ s"') 



Fig. 5. — Symmetric parameter SP as a function of the diffusivity rjs and the surface depth 
ds- The superposed hues indicate the contours of the dynamo cycle period over periods of 
years. Panel (a) shows the results for the slow meridional flow case {uq = 1000 cm s~^). 
Panel (b) shows the result for the fast meridional flow case {uq = 2000 cm s^^). 



